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20.  Abstract 


J.  Galambos  (1973)  for  modal  s v/i  th  '■bis  ornoertv.  Force,  for  most  functions 


F(x)  occuring  in  practice.  (7 


'b  Has  an  asymototic  distribution  H(x) 


as  n— > + oo.  This  asymptotic  distribution  is  necessarily  one  of  the  three 

classical  types  obtained  in  the  '920s  by  Fisher,  Tippett  and  Frechev. Hov/ever  , 

if  F (x)  is  not  known  (which  is  the  case  of  most  practical  situations),  then 

neither  the  normalizing  constants  a and  b , nor  the  actual  type  of  the  limit 

n n 

distribution  H(x)  can  be  determined.  Since  goodness  of  fit  tests  are  not 
sensitive  for  the  tail  of  the  distribution  F(x),  and  since  an  estimate  for 
F(x)  is  not  suitable  in  extreme  value  theory  (see  Galambos  (1978),  Example 
2.6.3).  some  other  techniques  are  needed  which  are  free  from  F (x) . The 
paper  is  devoted  to  such  a technique,  which  is  an  extension  of  the  approach 
by  J.  Pickands  III  (1975).  In  such  a method,  one  uses  a number  of  upper 
extremes  and  the  parametric  form  of  the  three  classical  types  for  H(x), 
known  as  the  von  Mises  form. 
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Let  X^,  X,,,...,  XR  be  random  variables  with  common 
distribution  function  F(x).  liat  = max(  X^  , X^ , . . . ,Xn ) . Assume 

that  the  interdependence  of  the  X..  is  such  that,  with  suitable 
constants  an  and  bn>  0,  P(Zn<an  + bnx)  is  asymptotically 
Fn(an  + b^x);  fioo  Chapter  3 in  the  book  by  the  present  author, 

J.  Gnlambos  (1978)  for  models  with  this  property.  Hence,  for 
most  functions  F(x)  occuring  in  practice,  (ZR  - a^/b^  has 
an  asymptotic  distribution  H(x)  as  n— >+oo.  This  asymptotic 
distribution  is  necessarily  one  of  the  three  classical  types 
obtained  in  the  1920s  by  Fisher,  Tippett  and  Frechet.  However, 
if  F(x)  is  not  known  (which  is  the  case  of  most  practical 
situations),  then  neither  the  normalizing  constants  and  b , 

nor  the  actual  type  of  the  limit  distribution  H(x)  can  be 
determined.  Since  goodness  of  fit  tests  are  not  sensitive  for 
the  tail  of  the  distribution  F(x),  and  since  an  estimate  for 
F(x)  is  not  suitable  in  extreme  value  theory  (see  Galambos  (1978) 
Example  2.6. 3) , some  other  techniques  are  needed  which  are  free 
from  F(x).  The  paper  is  devoted  to  such  a techninue,  which  is 
an  extension  of  the  approach  by  J.  Pickands  111  (1975)*  In  such 
a method,  one  uses  a number  of  upper  extremes  and  the  parametric 
form  of  the  three  classical  types  for  H(x),  known  as  the  von  Mine 
form. 
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IS 


A new  statistics!  method 


Let  Xlf  X2,...,  X^  be  identically  distributed  random 
variables.  Put  Zn  = max(X^,X2 » • • • »Xn) . Assume  that  there  are 
constants  an  and  bn>0  such  that  P(Zn<  an  + bnx)  can  be 
approximated  by  {^(a  + bnx),  where  F is  usually  the  common 

distribution  function  of  the  X’s.  This  implies  that  if  F has 
a smooth  enough  tail,  then  (Zn  - an)/bn  has  an  asymptotic 
distribution  II(x).  Our  aim  is  to  determine  the  form  of  H(x) 
in  the  case  when  F is  not  known.  The  approximation  of  F(x) 
by  the  empirical  distribution  function,  or  by  other  classical 
methods,  is  inefficient  in  extreme  value  theory;  see  Galambos 
(1978a)  for  detailed  discussion.  Hence  a direct  method  as  required 
to  get  information  on  H(x)  without  trying  to  determine  F. 

There  is  only  one  work  on  this  line  in  the  literature  by  J. 

Pickands  III  (1975).  Here  we  consider  an  alternate  model.  While 
Pickands  developes  a method,  in  which  data  will  automatically 
separate  the  form  of  H(x)  out  of  the  three  possibilities,  our 
method  is  a hypothesis  test,  in  which  we  make  a decision  on  the 
mentioned  form.  It  remains  an  unsolved  problem  to  compare  the 
two  methods  on  a theoretical  basis.  Because  the  nature  of  the 
two  approaches  is  different,  a clear  comparison  is  difficult. 
However,  advantages  and  disadvantages  of  one  method  or  the  other 
can  be  clarified  from  several  points  of  vi ew.  As  a matter  of 
opinion,  we  can  state  that  the  method  by  Pickands  would  be 
preferable  to  a test  if  the  accuracy  of  his  method  were  available. 
Since  the  Pickands  method  involves  three  different  passages  to 
limit  for  supporting  its  conclusions,  an  estimate  of  accuracy 
is  quite  important.  On  the  other  hand,  in  some  cases,  one  would 
need  the  actual  form  of  H(x),  for  which  a test  is  suitable, 
while  the  method  of  Pickands  leaves  this  undecided. 

The  method  described  below  is  in  an  early  stage  of  investiga- 
tion. However,  initial  results  and  some  numerical  computations 
suggest  that  it  is  a useful  attempt  and  it  thus  deserves  a thorough 
investigation. 


Basic  assumptions:  1./  The  X.  are  identically  distributed  with 

common  distribution  function  F(x);  2./  The  dependence  of  the 

X . is  such  that  the  distribution  of  2 can  be  approximated 

by  the  n-th  power  of  F(x)  for  large  values  of  3./  F(x)<l 

for  all  x;  4./  E(X.)  is  finite;  and  5./  There  are  constants  a 
» J n 

and  such  that  '^n""an^^n  ^as  an  ar,yraP^otic  distribu- 

tion H(x). 


The  above  assumptions  imply  that  there  are  constants  A 
and  B > 0 such  that  H(A+Bx)  equals  either  Q(x)  = exp(-e“x) 
or  c(x)  = exp(-x”c),  x>0,  c >0.  Our  aim  is  to  decide  by 
a statistical  method  based  on  the  X.  that  which  of  the  two  forms 
mentioned  applies. 


A test  for  the  form  of  H(x). 


Notice  that  the  X . 

0 


provide  a 


single  observation  on  Zn.Yet  a number  of  extremes  can  be  used 
to  provide  information  on  H(x).  The  actual  method  is  as  follows. 


Choose  a number 

»uch  X1,  s to 
J 


Select  those  X.  that  exceed  t.  Transform 
J 


Yi  * <Xj  - U/K(Xrt|x3>t)  . 


Test  that  the  distribution  of  these  Y.  is  unit  exponential 

J 

against  the  comnosite  hypothesis  that  this  distribution  belongs 
to  the  Pareto  family.  Accepting  the  hypothesis  means  acceptance 
of  H(A+Bx)  = H.  q(x),  while  rejection  of  the  null  hypothesis 

leads  to  the  acceptance  of  H,  (x)  as  the  limiting  distribution 

» c 

of  the  properly  normalized  maximum.  The  nower  of  the  test  depends 
on  the  choice  of  t,  which  choice  is  yet  to  be  investigated  in 
terms  of  some  optimality  criterion. 


The  theoretical  justification  of  the  above  procedure  is 
a classical  criterion,  combined  with  a result  of  L.  de  Haan  (1cj70), 
for  distributions  belonging  to  given  domains  of  attraction  for 
the  maxima.  Details  will  be  given  in  an  expanded  form  of  this 
present  paper. 
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